Ground state of the frustrated Hubbard model within DM FT: 
energetics of Mott insulator and metal from ePT and QMC 

N. Blumer a '* E. Kalinowski b 

a Institut fur Physik, Johannes-Gutenberg-Universitat, 55099 Mainz, Germany 
b FB Physik, Philipps Universitat, 35032 Marburg, Germany 



Abstract 

We present a new method, ePT, for extrapolating few known coefficients of a perturbative expansion. Controlled by 
comparisons with numerically exact quantum Monte Carlo (QMC) results, f th order strong-coupling perturbation 
theory (PT) for the Hubbard model on the Bethe lattice is reliably extrapolated to infinite order. Within dynamical 
mean-field theory (DMFT), we obtain continuous estimates of energy E and double occupancy D with unprece- 
dented precision 0(10 ) for the Mott insulator above its stability edge U c i ~ 4.78 as well as critical exponents. 
In addition, we derive corresponding precise estimates for E and D in the metallic ground state from extensive 
low-temperature QMC simulations using a fit to weak-coupling PT while enforcing thermodynamic consistency. 
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The half-filled Hubbard model has been subject of 
intense research; after clarifying the magnetic phases, 
DMFT studies have lately focused on the paramagnetic 
Mott metal-insulator transition observed in frustrated 
versions. In the following, we shortly summarize recent 
work; for more complete and self-contained presenta- 
tions, see Ref. [1] and upcoming publications. 

Kato-Takahashi perturbation theory (PT) for the 
half-filled frustrated Hubbard model 
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"Bethe" density of states (W = At) yields [1] the energy 
E{U) = E PT (U) + 0(t 12 I!' 11 ), where for t = 1: 
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Both the 10 th order (in t) expression (1) for £pt and 
the resulting estimate Dpt = dEpT / dll for the double 
occupancy are converg ed within 0(1O" 5 ) for U > 6. 
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Fig. 1. Energy E and double occupancy D of Mott insulator: 
PT of various order (lines) in comparison with QMC (circles). 

Similar precision can be reached at finite temperatures 
T using an advanced quantum Monte-Carlo scheme [1]; 
due to the almost negligible temperature dependence 
in the insulating phase this also applies to the Mott in- 
sulating ground state. As seen in Fig. 1, both methods 
indeed agree very well for U — 6 (and deviate below). 

The idea of the extended perturbation theory (ePT) 
is to extrapolate coefficients in the PT series £W = 
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Fig. 2. Construction of ePT: Extrapolation of PT values for 
t/ c [n] (squares) defines coefficients for smaller 1/n (circles). 
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Fig. 3. Mott insulator: Energy/double occupancy from QMC 
(circles), ePT (solid line), and truncated ePT (dashed/dotted 
lines); all results after subtracting 10 th order PT. 

X/fca a 2i U x ~ 2 ' by fitting coefficient ratios f/ c i[2i] = 
\J a,2i+2/a2i to U c i[n] U c i + ui rT 1 +112 rT 2 . As seen 
in Fig. 2, the convergence to this assumed asymptotic 
form is fast; consequently, both the extrapolation to 
Ud = limj_>oo t/ c i[2i] and the extraction of higher- 
order coefficients f7 c i[2i] are reliable. We can further 
derive algebraic exponents: a n oc n T U"i, E CT i t (U) oc 
{U - U cl ) T ~\ D CIlt (U) oc (17 - U cl ) T ~ 2 . Analytic 
and numerical arguments support [1] r = 3.5, U c i = 
4.78. Fig. 3 shows the perfect agreement between ePT 
and QMC. These results mark a breakthrough in ac- 
curacy for correlated electron systems; the ePT curves 
(for which parameterizations are given in [1]) should 
be/are used as authoritative benchmarks [2]. New ev- 
idence confirms our estimate of U c i = 4.78 which in 
turn strongly supports our analysis of critical expo- 
nents (for which no reliable comparisons exist). 

Continuous ground state estimates are much more 
difficult to obtain from finite-temperature methods 
(such as QMC) in the metallic phase due to the much 
stronger temperature dependencies (with divergent 
linear coefficient of specific heat for U — * U C 2 ~ 5.84). 
On the basis of careful extrapolations of data for 
l/T = 15, 20, 25, 32, 40 (Fig. 4) we have obtained such 
estimates with precision O(10 -4 ) as shown in Fig. 5; 
this accuracy suffices in order to discriminate between 
QMC estimates (circles) and 2 nd order weak-coupling 
PT (dotted lines) . The continuous and thermodynam- 
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Fig. 4. Metal: QMC Energy/double occupancy (symbols) vs. 
squared temperature T and extrapolation T — > (lines). 
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Fig. 5. Ground state of frustrated Hubbard model within 
DMFT: Energy E and double occupancy D of metal and in- 
sulator match at critical point U c 2 ~ 5.84 (circle). 2 nd order 
weak-coupling PT (dotted lines) shown for comparison. 

ically consistent fits (dashed lines) to QMC in Fig. 5 
include the well-known 2 nd order term and fits to few 
higher even order coefficients under the constraint of 
matching the ePT estimates for the insulator at U = 
U C 2 ~ 5.84 (circles). Note that our implicit estimate 
of the (weak-coupling) 4 th order coefficient deviates 
strongly (and has opposite sign of) the analytic result 
found in the literature [3]. Explicit parameterizations 
are deferred to a full-length publication which will also 
discuss additional observables implicit in our deriva- 
tions (such as the linear coefficient of the specific heat) . 
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